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We calculate angle-dependent magnetoresistance oscillations (AMRO) for interlayer transport of 
cuprate superconductors in the presence of (n, i r) order. The order reconstructs the Fermi surface, 
creating magnetic breakdown junctions; we show how such magnetic breakdown effects can be incor¬ 
porated into calculations of interlayer conductivity for this system. We successfully fit experimental 
data with our model, and these fits suggest a connection between (7r, 7r) order and the anisotropic 
scattering observed in overdoped cuprates. This work paves the way for the use of AMRO as a tool 
to distinguish different kinds of ordered states. 


I. INTRODUCTION 

Understanding the nature of broken symmetry phases 
in the thermodynamic phase diagram of the cuprates 
is a key step toward understanding the origin of high- 
temperature superconductivity. For example, the dis¬ 
covery of the pseudogafpl has fueled the search for many 
kinds of order®!] including nematic phases that could 
strongly enhance tJ®. Yet of the broken symmetries 
connected to unconventional superconductivity, antifer¬ 
romagnetism remains one of the most important, appear¬ 
ing in cuprate, iron-pnictide, organic, and heavy fermion 
materialise. 

Unambiguous evidence of the presence of Fermi sur¬ 
face reconstruction arising from broken symmetry or¬ 
der has come from quantum oscillation measurements in 
both hole-doped 13 AS and electron-doped-SO cuprates at 
low temperatures and high magnetic fields, but the na¬ 
ture of the broken symmetry remains a matter of consid¬ 
erable debate. Antiferromagnetic (7r, 7t) reconstruction 
has been prop osed for both the hole- and electron-doped 
material :] 16 f l7 ( but recent evidence for a (possibly field- 
induced) charge density waviSIK 21 h as suggested more 
complex orders are driving the reconstruction. 

The ability to experimentally differentiate between 
these different ordered states is crucial. In this work, 
we suggest that interlayer angle-dependent magnetoresis¬ 
tance oscillations (AMRO) can be used to distinguish dif¬ 
ferent kinds of long-range ordered states in the cuprates. 
Angle-dependent magnetoresistance is a sensitive probe 
of the Fermi surface of a materiaP^HMI anc j can therefore 
be used to investigate the geometry of a reconstructed 
Fermi surface. The measurement is also sensitive to the 
energy scale of any (translational) symmetry-breaking 
ordei 27 . This energy scale is related to a “magnetic 
breakdown field,” as we will describe below. Importantly, 
these effects on AMRO can be observed even in materi¬ 
als that do not show quantum oscillations. Thus, the 
measurement is useful in systems in which sample disor¬ 
der is high , or in which the order has a small correlation 

lengtlPAO. 

AMRO data from TUE^CuOe+a provided the ear¬ 


liest transport evidence for the existence of a three- 
dimensional Fermi surface in an overdoped cuprat^^. 
The temperature evolution of the AMRO is consistent 
with a superposition of isotropic and anisotropic scat¬ 
tering rates about the Fermi surfac^®, and it was de¬ 
termined that these do not have the same temperature 
dependence: the isotropic scattering rate is quadratic 
with temperature (as expected of an ordinary Fermi liq¬ 
uid) while the anisotropic scattering rate is linear (con¬ 
necting it to the non-Fermi liquid physics of the cuprate 
phase diagram). Additionally, the anisotropic scattering 
is strongest in the anti-nodal region of the Fermi sur¬ 
face. Therefore, it has been suggested that the anoma¬ 
lous scattering temperature dependence may be related 
to (7r,7r) fluctuation^®®} possibly originating from an- 
tiferromagnetism. We show that the effects of such or¬ 
der on the AMRO can indeed make a natural connection 
with the physics of the observed scattering anisotropy, 
even though there is no static order in the system. 

We demonstrate this connection by simulating the 
AMRO of a model cuprate material in the presence of 
antiferromagnetic order. The interpretation of AMRO 
measurements requires efficient and versatile calculations 
of the magnetotransport of a given Fermi surface so that 
models can be compared to experimental results. These 
calculations are more challenging in the presence of static 
order that reconstructs the Fermi surface. We have devel¬ 
oped a general method to perform such calculations for 
quasi-two-dimensional (Q2D) materials, based on previ¬ 
ous work in organic metals®. It is both easy to imple¬ 
ment and computationally inexpensive. In Section |TT| we 
use this method to calculate the interlayer magnetore¬ 
sistance of a tetragonal Q2D material in the presence of 
(7r, 7 r) order and including the effects of magnetic break¬ 
down. In Section |III| we apply this model to the known 
Fermi surface of TUf^CuOe-M^and show that the tem¬ 
perature dependence of the AMRO can be captured by 
this magnetic breakdown model. In Section |TV| we discuss 
the physical consequences of this model and its potential 
range of applicability for distinguishing between differ¬ 
ent kinds of broken symmetry order in the cuprates. Our 
general method is laid out in detail in Appendix [A] 
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II. AMRO IN THE PRESENCE OF ( 77 , 77 ) ORDER 

As a first application of our method, we wish to un¬ 
derstand how the AMRO of the cuprates is affected by 
antiferromagnetism. We therefore consider the case of 
a Q2D tetragonal material under static (7r, 7 r) antiferro¬ 
magnetic order (though the model below can be applied 
to any ( 77 ,7r) order). As shown in Figurejlja), the original 
Brillouin zone of such a material will have a square cross- 
section with primitive reciprocal lattice vectors along k x 
and k y ; we define all azimuthal angles in this paper with 
respect to k x . In the presence of ( 77 , 77 ) antiferromagnetic 
order, the Brillouin zone is halved in cross-section, re¬ 
sulting in a reconstruction of the Fermi surface as shown 
in Figure |T](b,c) . This new reconstructed Brillouin zone 
will have primitive reciprocal lattice vectors along k' x and 
ky, which are rotated by 45 °with respect to k x and k y . 



FIG. 1. (Color online) Fermi surface reconstruction of a Q2D 
material under ( 77 , 77 ) antiferromagnetic order, as viewed along 
k z . (a) The Fermi surface (FS) and first Brillouin zone (BZ) 
of a Q2D tetragonal material; (b) the reconstructed BZ and 
(c) reconstructed FS of the material following the onset of 
( 77 , 77 ) antiferromagnetic order; (d) the repeated-zone view of 
the reconstructed FS, illustrating the small cross-sections of 
FS that have replaced the unreconstructed cylindrical FS. The 
grey line illustrates Bragg diffraction between two magnetic 
breakdown junctions at the BZ boundary. The angle £ is also 
defined here; it will be used in our conductivity calculations. 

Quasiparticles traversing the Fermi surface will Bragg 
diffract at the reconstructed Brillouin zone boundaries, so 
they will travel along three distinct Fermi surface pock¬ 
ets as shown in Figure m- However, in a large mag¬ 
netic field, the quasiparticle path in real space may be 


curved sufficiently to avoid Bragg diffraction. This is 
known as magnetic breakdown (MB), and can be thought 
of as a tunneling in fc-space from one pocket to the 
nextPS. The probability to tunnel in this way is given 
by p = e ~ B °! B , where Bq is the breakdown field and is 
a material-dependent constant proportional to the gap 
in fc-space between Fermi surface sectional. At every 
instance the quasiparticle path reaches a Brillouin zone 
boundary, the quasiparticle may either Bragg diffract or 
undergo MB; thus, such points in the quasiparticle path 
are known as MB junctions. 

We must take the effect of these MB junctions into 
account when calculating conductivity. The conductivity 
of a Q2D material in a magnetic field can be calculated 
using the Shockley tube integral form of the Boltzmann 
transport equation^, 
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where ipo is the initial azimuthal position of the quasipar¬ 
ticle and if is its position after some time t has passed^. 
The effective mass of the quasiparticle is represented by 
m* and the cyclotron frequency is oj c = eB/m*. The 
velocities in Eq. [l] are Fermi velocities. 

The vector k B points parallel to the magnetic field and 
defines the orbital path of a quasiparticle. We integrate 
across all values of its magnitude. For a given magnitude, 
the tip of the vector will touch a single quasiparticle orbit 
which can be defined by /c?, the /^-position of the orbital 
plane at the center of the Fermi surface (see Figure [2j) . 
The magnetic field’s direction is defined by a polar an¬ 
gle (0) with respect to k z and an azimuthal angle (<j>) 
with respect to k x . We can write \k B \ = k®cos(6) and 
therefore convert our integral over k B to one over k®. 

Since AMRO is a probe of interlayer conductivity, we 
want to calculate o zz . This means we need an expression 
for v z , which we can obtain from a symmetry-constrained 
model of the Fermi surface. The following equation de¬ 
scribes a Q2D Fermi surface of a layered tetrag onal ma¬ 
terial with simple cosine warping along feJUMl 


E F (k z ,f) 


H 2 k^(f) 

2m* 



( 2 ) 


In the above, t± is the interlayer hopping, while kp(tp) 
and F(f) parameterize the Fermi surface in the az¬ 
imuthal cylindrical coordinate. The in-plane and out- 
of-plane lattice parameters are denoted by a and c, re¬ 
spectively. Using v z = £ , we find the interlayer 

velocity to be 

v z( k , <p) = sin F (<P)- (3) 
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FIG. 2. When a magnetic field is applied to a Q2D ma¬ 
terial, quasiparticles will trace out orbits on the cylindrical 
Fermi surface that are perpendicular to the applied field. On 
the lowest orbit in this figure, we define the parameters k° z 
and ks- On the upper two orbits, we illustrate the fact that 
a quasiparticle undergoing Bragg diffraction moves to a dif¬ 
ferent cross-section of the Fermi surface. The dashed lines 
represent the Brillouin zone boundaries for the reconstructed 
Fermi surface. Note that the azimuthal angle, <f>, is defined 
with respect to k x of the unreconstructed system. 


When a quasiparticle Bragg diffracts at the Brillouin 
zone boundary it will have a momentum change given 
by a reciprocal lattice vector, so its momentum in the 
^-direction will not change. The particle will jump to a 
different “slice” of the Fermi surfac^H 0 f k z (see Figure 
[2]), changing its value but preserving k z . The amount 
by which k z changes after Bragg diffraction depends on 
which MB junctions are involved; for each pair of MB 
junctions, the value of A k z can be calculated using purely 
geometric means (see AppendixjC]) . Therefore, we can see 
that for a given quasiparticle, 


k z (ip) = k z — kp(ip) tan(0) cos (ip — </>) + E rij(ip)Ak { z J \ 

3 

(4) 

where A kf'P are the possible changes of k z from Bragg 
diffraction and rij are the number of times each has 
occurecP^. Note that we have neglected the influence of 
t _l on particle motion, which is a reasonable ommission 


except for 6 « 90°. Setting rij(<po) = 0 we find 
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up to a constant of proportionality. 

Performing the integration over k z , we arrive at 
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where we define G(ip) = § • kp(ip) tan(0) cos (ip — <j>). 

We neglect the constant prefactor and use 

cos(a;) = Re[e“] to write 

/• 2 tt p 00 

I dipo F((p 0 ) / dip F(ip)x 

0 Jipo (7) 
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Note that the value of the integrand changes whenever 
the quasiparticle undergoes Bragg diffraction, due to the 
term n j(ip) A ki^. In order to evaluate the integral, 
we must be able to account for all possible trajectories of 
each quasiparticle. 

Following the method of Nowojewski et alW^ 1, we will 
separately consider the motion of quasiparticles starting 
in the 8 different segments of the Fermi surface, then 
sum their contributions to the conductivity. To do so, 
we rewrite the above integral in a vectorized form: 
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In this equation, the dot product with sums up 
all the possible initial positions of the quasiparticle, Ai n i t 
describes the initial motion of the quasiparticle up to 
an MB junction, and X v describes the contribution to 
conductivity when the particle is between MB junctions. 
Each vector is 8-dimensional, and they are defined as 
follows: 
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where g = 1 — p and we have defined a = e _2 ^/ aJcT and 
b = e _ ( 7r / 2_2 £)/‘‘ JcT . See Appendix | a| for an explanation 
of the elements of T. 

As a simplification, we have assumed that the gaps 
that open in the Fermi surface upon reconstruction are 
of a negligible length in k- space: we take the MB junction 
that ends one section of the Fermi surface to be in the 
same position as the MB junction that begins the next 
section. 

With a zz in this vectorized form, we can quickly cal¬ 
culate numerical values for the conductivity with varying 
9 and </>. 


III. APPLICATION TO A CUPRATE 
SUPERCONDUCTOR 

We are now in a position to apply this model 
to a real system. We focus on Tl 2 Ba 2 Cu06-t-,5, 
since this is the cu prate tha t has been studied 
the most with AMR.O ^ 1 29 1 36 1 37 1 . As described by 
Hussey et alW^ the Fermi surface of Tl 2 Ba 2 Cu06+,5 
can be parameterized by kp(tp) = fc 0 o + feio cos(4<p) and 
F(ip) = fc 2 i sin(2p) + kei sin(6<p) + fcioi sin(10</3). The 
coefficients k mn label an expansion of the Fermi surface 
in cylindrical harmonics appropriate for the space group 
symmetry of this materiaPS. 

The AMRO of an unreconstructed Fermi surface can 
be produced by setting Bq = 0, as shown in Figure 
[3] (a-d) for two convenient values of uj c t. Note that 
(jj c t = 0.5 reproduces the experimental AMRO observed 
by Hussey et al. at 4.2 K 2 S. The AMRO for a system 
with (7r,7r) antiferromagnetic order is shown in Figure [3] 
(e-h). This shows many qualitative differences with the 
unreconstructed state. The peak at 9 = 0 is strongly sup¬ 
pressed in the reconstructed Fermi surface. In addition, 
there are more Yamaji angles (peaks in the AMRO) for 


is a vector giving the azimuthal position of each MB junc¬ 
tion and the angle £ is defined in Figure [ljd). 

The matrix T accounts for the connections between 
orbit segments, as well as the exponential damping of 
the integrand upon traversing a segment of Fermi surface. 
For our system, it is an 8 x 8 matrix: 
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low polar angles 9 in the unreconstructed state than the 
reconstructed state. 

The evolution of the AMRO as we go from b o/b ss 0 to 
b o/b = 10 for oj c t = 0.5 bears a striking resemblance to 
the evolution of the AMRO in Tl 2 Ba 2 CuOg+a with in¬ 
creasing temperat ure, m ost notably the disappearance of 
the hump at 9 = (pB 39 1 This seems surprising given that 
Tl 2 Ba 2 Cu 06+5 is not known to exhibit any static anti¬ 
ferromagnetic order, though it has been shown to have 
strong antiferromagnetic fluctuation^^]. 

We explore the possibility that our AMRO calcula¬ 
tions can capture some of the physics of Tl 2 Ba 2 Cu06+<5. 
Using the form of er 2Z above, we have produced simu¬ 
lations of out-of-plane resistivity as a function of angle 
using existing data for a sample with T c = 15 K reported 
in Ref. [57[ The low-temperature (4.2 K) AMRO of 
Tl 2 Ba 2 Cu 06+5 is well-fit by a simple model with no an¬ 
tiferromagnetic order (Bq = 0), and using this data the 
functions fc|,(<p) and F(ip) that describe the Fermi surface 
can be fully determined in good agreement with previous 
worl J 28 | 36 lPp See Appendix § for more information on 
our determination of these parameters. We used geomet¬ 
ric methods to solve for £ and A ki^ in this system, as 
explained in Appendices [B] and [C] respectively. To study 
the temperature-dependent AMRO above 4.2 K we al¬ 
lowed only two free parameters: w c r and Bq. Note that 
in contrast to Ref. |37l w c r is fixed to be isotropic with 
azimuthal angle ip. We ran simulations across a large 
range of parameter space and used a least-squares fitting 
approach to determine the values of w c r and Bq at each 
temperature. Our best fit to the data is shown in Fig¬ 
ure |4j showing excellent quantitative agreement with the 
AMRO of Tl 2 Ba 2 Cu06 + a. The temperature dependence 
of w c r and Bq can be extracted from these simulations, 
and these are shown in Figure [5] 

As can be seen in Figure [3j more features are appar¬ 
ent in the AMRO when w c r is higher, making it eas- 
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FIG. 3. (Color online) Calculated dimensionless resistivity of Tl2Ba2Cu06+5 as a function of the orientation of applied magnetic 
field for different values of b o/b and lo c t. The resistivity in plots (a,b) and (e,f) is scaled up by a factor of 5 relative to the 
remaining plots, for visual clarity. (a,c,e,g) Interlayer resistivity as a function of both 6 and rj>. (b,d,f,h) Interlayer resistivity 
as a function of 9 for select values of (j ); these azimuthal angles were chosen to facilitate comparison to Ref. 1251 


ier to distinguish the effects of changing B 0 . If oj c t is 
decreased (by lowering magnetic fields, raising tempera¬ 
tures, or lowering sample quality), each quasiparticle will 
traverse less of the Fermi surface before it scatters. As 
illustrated in Figure 6 of Ref. HU this causes the ampli¬ 
tude of AMRO to be reduced, which makes an accurate 
determination of Bq more difficult. Thus, the error of 
our fitting parameters is greater at higher temperatures. 
Indeed, Ref. includes AMRO data taken at 90 K and 
110 K, but we were not able to accurately determine Bq 
at those elevated temperatures. 


While Ref. [37] reproduced the observed AMRO using 
an anisotropic scattering rate, we find good quantita¬ 
tive agreement with the data using a magnetic break¬ 
down model with an isotropic scattering rate. Bragg 
scattering at a MB junction mimics the effect of an 
anisotropic scattering rate on an unreconstructed Fermi 
surface. However, importantly the magnetic breakdown 
model connects specific parts of the Fermi surface in 
a single (Bragg) scattering event, while the anisotropic 
scattering rate is a broad modulation of the quasiparti¬ 
cle lifetime about the Fermi surface. The similarity of 
the two models in reproducing the AMRO suggests that 
the apparent anisotropic scattering rate is a symptom of 
antiferromagnetism, perhaps involving short-range fluc¬ 
tuations. This could explain the different temperature 
dependence of the isotropic and anisotropic components 
of the scattering rate observed in Ref. [2SJ 


IV. DISCUSSION 

The behavior of lo c t in Figure [5](a) indicates a nat¬ 
ural (approximately linear) increase in the scattering 
rate with temperature. The evolution of Bq may reflect 
deeper physics. As shown in Figure [5jb) , B 0 increases 
quickly with temperature, peaking around 45 K. The pa¬ 
rameter Bq is a measure of the probability of Bragg scat¬ 
tering. For a static reconstructed Fermi surface, this is 
related to the separation between reconstructed sections, 
which is in turn proportional to the bandgap 27 . There¬ 
fore, under static reconstruction we would expect Bq to 
be largest at 0 K and decrease weakly with increasing 
temperatur d 27 * 42 ! In the presence of antiferromagnetic 
fluctuations, similar scattering events might still occur 
at points where the reconstructed Brillouin zone inter¬ 
sects the Fermi surface. In this case, B 0 will play two 
roles: in addition to parameterizing the separation be¬ 
tween sections of Fermi surface, it also reflects the prob¬ 
ability of Bragg scattering within the time/length-scale of 
the fluctuation^^. Note that in the overdoped cuprates, 
there is a known crossover in the transport from Fermi 
liquid- to non Fermi liquid-like behavior with increasing 
temperature that is thought to be associated with critical 
fluctuation^^. j n this picture, the increase of B 0 with 
temperature (Figure |5j[b)) can be interpreted as an in¬ 
crease in antiferromagnetic fluctuations. As the tempera¬ 
ture rises and antiferromagnetic fluctuations grow, quasi¬ 
particles have a non-zero chance of undergoing Bragg 
diffraction when they reach MB junctions, so Bq attains a 
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For the above to be plausible, the antiferromagnetic 
fluctuations of the system should be on a long enough 
timescale to affect the quasiparticles’ motion about the 
Fermi surface: the timescale of an antiferromagnetic 
fluctuation should be longer than the time it takes for 
a quasiparticle to traverse a section of Fermi surface 
from one MB junction to the next. The antiferromag¬ 
netic fluctuations in I^-xSr^CuC^ near optimal doping 
have a frequency that is roughly linearly proportional to 
temperature^. Taking this as a guide, we estimate that 
the timescale of an antiferromagnetic fluctuation will be 
of the order taf ~ ttt- Meanwhile, the time for a quasi¬ 
particle to cross the smallest section of Fermi surface be¬ 
tween two MB junctions is given by tqp ~ ■ 

Therefore, our condition taf > tqp is equivalent to 


FIG. 4. (Color online) Temperature dependence of the in¬ 
terlayer AMRO of overdoped Tl2Ba2Cu06+s (T c = 15 K) 
at a fixed field of 45 T and a fixed azimuthal direction of 
0 = 7°. The solid lines are c-axis magnetoresistivity data, 
taken from French et al These data have been normalized 
to the zero-field resistivity of the sample at each temperature. 
The dashed lines are simulations of AMRO for Tl2Ba2Cu06+a 
under antiferromagnetic order, calculated as described in the 
text. When producing these simulations, all the parameters 
related to Fermi surface geometry were fixed and the only 
parameters allowed to vary with temperature were uj c t and 
B 0 . 
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FIG. 5. (Color online) Temperature dependence of a) x /u,,t 
and b) Bo extracted from fits to data. Error bars are stan¬ 
dard errors extracted from the covariance matrix of the least- 
squares fitting at each temperature (see Appendix 0- 


Hue > • k B T. (12) 

For this system the requirement is approximately 

^>3.6x10 9 s~ 1 K~ 1 . (13) 

Using m* « 5 we find w c « 1.6 x 10 12 s -1 at 45 
T. Therefore, antiferromagnetic fluctuations could be ex¬ 
pected to affect quasiparticle motion up to T ~ 400 K, 
much higher than the temperature regime studied in this 
paper. 

The magnetic breakdown picture of the effect of an¬ 
tiferromagnetic fluctuations on AMRO could be sub¬ 
stantially improved by including a more realistic model 
of the MB junctions in a fluctuating system that in¬ 
cludes, for example, a distribution of ordering wavevec- 
tors about (n, 7r) 17 . Nevertheless, this simple model cap¬ 
tures many of the important features observed in the 
temperature-dependent AMRO without the need for a 
multi-component scattering with a differe nt nodal and 
anti-nodal temperature dependenc d 29 | 36 | 37 J . Indeed, the 
success of this simple model indirectly connects this un¬ 
usual scattering anisotropy to antiferromagnetism, which 
may be useful for understanding other transport prop¬ 
erties (see Appendix E]). Moreover, our results suggest 
there is a potential link between Bq and the dynamic sus¬ 
ceptibility Imy. If this connection can find a sound the¬ 
oretical basis, it may open the way for the use of AMRO 
as an experimental probe of magnetic scattering. 


non-zero value. At still higher temperatures, the antifer¬ 
romagnetic correlation time is so short that the effect of 
Bragg scattering decreases, resulting in a decrease in Bq. 
The evolution of Bq looks strikingly similar to the evolu¬ 
tion of the imaginary part of the dynamic susceptibility 
Imy; (which is a measure of the magnetic scattering) seen 
in a number of neutron experiments in cuprate supercon¬ 
ductors; consider, for example, Figure 10 of Ref. 05j We 
therefore suggest that the temperature dependence of B 0 
in Figure[5](b) reflects the effect of antiferromagnetic fluc¬ 
tuations on the magnetotransport. 


V. CONCLUSION 

We have developed a simple and computationally in¬ 
expensive numerical method to calculate AMRO in lay¬ 
ered two-dimensional materials with ( 7 r, 7 r) antiferro¬ 
magnetic order. This model can be applied to both 
hole- and electron-doped cuprates with an appropri¬ 
ately adjusted Fermi surface parameterization for di¬ 
rect comparison with experimental data. In addition, 
our numerical method can easily be applied to or¬ 
dered states other than antiferromagnetism, such as the 
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charge-ordered states re cently proposed in underdoped 
YBa 2 Cu 3 0 6 + 1 ^ 21 — ^ and HgBa 2 Cu0 6+( P. We have 
shown that an antiferromagnetic Fermi surface recon¬ 
struction with a temperature-dependent magnetic break¬ 
down field can fit the AMRO of Tl 2 Ba 2 Cu06+a, an 
overdoped compound with no static order. The agree¬ 
ment between our fits and the AMRO data suggest 
that th e appare nt scattering anisotropy observed in these 
system j ^ 9 * 36 ! 37 ! is connected to antiferromagnetic fluctua¬ 
tions, and indeed that the MB field, can potentially 
be used as an experimental measure of such fluctuations. 
This would make AMRO a good complement to scatter¬ 
ing probes of fluctuations, such as neutron scattering and 
resonant inelastic X-ray scattering. We propose that fu¬ 
ture AMRO experiments at higher magnetic fields and in 
materials where Imy has been determined independently 
by neutron scattering would provide an instructive com¬ 
parison to test the validity of this connection. 
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Appendix A: General method for calculating 
conductivity with magnetic breakdown 

In this Appendix we describe a step-by-step method 
to calculate a zz in a Q2D material with magnetic break¬ 
down effects. 

1. Consider the full, warped-cylindrical Fermi surface 
that would exist were the Fermi surface not recon¬ 
structed. Using existing data or theories, determine 
a likely form of this Fermi surface as a function of 
k z and p. This may be exactly fixed or it may 
contain free parameters to be fitted. 

2. Use the Fermi surface to determine v a (k,p) and 
vp(k, p) for the element cr a p in question. Note that 
v x and v y are not simply proportional to k x and k y 


for a noncircular Fermi surface; see the section on 
in-plane transport below. 

3. Insert these velocities into the Boltzmann transport 
equation as given in Equation 1 of the main text. 
Wherever k z appears in the integrand, replace it 
with the following function of tp: 

k z (p) = k z — k\,{p) tan(0) cos(y> — </>) + 5>i {p)Ak^ 

(Al) 

4. Replace the integral over with an integral over 
k z multiplied by cos (9), and perform the integra¬ 
tion over kP z . At this point, it should be possible 
to write the Boltzmann transport equation in the 
form 

/»27T pOO 

a ap = C x / dp 0 Mpo) / dp f 2 {p)e c2 ^ "'<*>**« 

JO J tp o 

(A2) 

where /i and / 2 are functions, and Ci and C 2 are 
constants. Note that C 2 will be zero if /? = x or y. 

5. Determine geometrically where the Fermi surface 
will intersect the (reconstructed) Brillouin zone. 
These points are the magnetic breakdown junc¬ 
tions. Write a vector, AT, giving the azimuthal 
position of each junction and ending at the loca¬ 
tion of the first junction plus 27 t. Be sure that the 
definition of p = 0 for this vector is consistent with 
the definition of p = 0 for the Fermi surface warp¬ 
ing. The length of AT will be n + 1, where n is the 
number of MB junctions around the Fermi surface. 

6. Define three vectors of length n as follows: 

fMj+l 

A V0 \j] = e ~ M ^^ dpofi(po) 

JMj 

, r M i +1 

A v \j] = e^cT / h{(p) (A3) 

r M j +1 

Ainit [j] = / d( , f 

J vo 

7. Define the n x n matrix T. Each row (column) 
of r corresponds to a specific section of the Fermi 
surface between two MB junctions. The first row of 
T corresponds to the section between the first and 
second MB junctions, as defined in the vector AT; 
the second row corresponds to the section between 
the second and third MB junctions; and so on. The 
elements in each row are as follows: 






0 


if section i has no connection to section j 


dip 


r ij — ' 


if section i is connected to section j through magnetic breakdown 

(i.e., the quasiparticle on section i goes to section j by following the 

full cylindrical Fermi surface) (A4) 


cij( 1 — p)e C2Ak * +1) ^ J if section i is connected to section j through Bragg diffraction 

(i.e., the quasiparticle on section i goes to section j by following the 
reconstructed Fermi surface) 


where p = e~ B °/ Bcos \ e ) is the magnetic breakdown 
probability, and a,i = The term 

cii accounts for the damping of our integrand as the 
quasiparticle traverses the i th section of the Fermi 
surface. Note the term A : after traversing 
the i th section of the Fermi surface, the quasiparti¬ 
cle would Bragg diffract from the (i + l) th magnetic 

breakdown junction. The terms A/ci ^+1 ' l ~ i ’" , can be 
calculated as described in the section “Calculating 
A ki j)v below. 

8 . Using the objects defined above, calculate the con¬ 
ductivity for a given direction of the applied field: 

<Tap{e,<l>) = C 1 - Re[A„ 0 • (Ainit + IXI-r)-^)] (A5) 

Note that the dot product A Vo • Aj r ,i t yields a double 
integral over tpo and ip and must be evaluated as 
such. 



FIG. 6. (Color online) The angle £ is defined with respect 
to the reconstructed Fermi surface of Tl2Ba2Cu06+<5. The 
Fermi surface is shown in the repeated zone scheme, with the 
reconstructed Brillouin zone overlaid. 


Appendix B: Calculating £ 


Appendix C: Calculating A 


The angle £ is defined as shown in Figure [6j If the 
Fermi surface were completely cylindrical, it would obey 


c °s(£) 


7r 

oaf kF 


where oaf is the in-plane lattice parameter of the an- 
tiferromagnetically ordered system and kp is the Fermi 
momentum. We may neglect the interlayer warping of 
the Fermi surface, which is relatively weak, but not the 
in-plane warping. Therefore, we have the relation 


cos(£) = 


7r 

aAF(k 00 + k 4 o cos(4£)) 


which can be solved self-consistently for £. We know 
that oaf = t/2 a. We use a = 0.3866 nm, as given 
by Analytis et alW- 1. We use /coo = 7-30 nm -1 and 
/c 4 o = —0.234 nm -1 . These are the values found by 
French et al. from fitting their 4.2 K AMRO dataP^ 
and they are consistent with the results of our fits (see 
above). Using these values, we find £ « 40.18°. Due to 
uncertainty in the Fermi surface fits, we cannot calcu¬ 
late £ with accuracy beyond two significant digits. We 
therefore round to £ = 40° for use in our fits to high- 
temperature data. 


As stated in the main text, we can define a vector giv¬ 
ing the azimuthal position of each magnetic breakdown 
(MB) junction as follows: 

nr nr nr »7T »7T 

M=--£+[0, 2£, -, —+2£, 7r, 7 t+ 2£, y +2£, 2 tt] 

^ (Cl) 

The position of these MB junctions on the (unrecon¬ 
structed) Fermi surface is shown in Figure [T] 

To find the values of Ak9\ we must know where a 
quasiparticle goes when it experiences Bragg diffraction 
at a given MB junction. To determine this, we need 
only see which MB junctions are connected by reciprocal 
lattice vectors of the reconstructed Brillouin zone. They 
are the following: 1 6, 2 5, 3 8, 4 f> 7. 

An easy way to understand these pairings is to examine 
the small Fermi surface orbits that the quasiparticle will 
follows if it Bragg diffracts at every junction (see Figure 

©• 

As stated in the main text, a quasiparticle undergoing 
Bragg diffraction in this system will have ki^ = ki^. 

We neglect the weak interlayer warping of the system; 
then for a quasiparticle on a particular slice of the Fermi 
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FIG. 7. (Color online) Position of the 8 magnetic breakdown 
junctions on the Fermi surface of Tl2Ba2Cu06+« under (tv, n) 
order. The reciprocal space axes shown correspond to the 
reconstructed Brillouin zone. 

surface, we can write 

k z (p) = k° z — kp(ip) tan(0) cos (p — </>). (C2) 

This leads to the condition 

tan(0) cos (<fi -</>) = 
k^' 1 — kp((pf) tan(0) cos (ipf — (j>) 

and therefore 

A £■(*—*•/) — h°(/) _ 1,0(i) _ 

2, - z ^Z - 

tan(6)[kp((p f ) cos(<pf—(j)) - fc|(^)cos(^ - </>)]. 

(C4) 

Since tpi and ipf are given by the vector M , we now have 
everything we need to solve for for each possible 

Bragg diffraction. For example, if a quasiparticle is going 
from MB junction 1 to MB junction 6 we have ipi = f — £ 
and Pf = \ — £ + 7r + 2t; = ^+(;. We can use these to 

solve for A/ci 1_>6 \ which we denote as Afci 1 ' 1 for the sake 
of brevity. 


Simultaneously with fitting the Fermi surface geome¬ 
try, we used the 4.2 K data to fit for the misalignment 
of the crystal with respect to the magnetic field; see 
Analytis et al. for details on the significance and cal¬ 
culation of this misalignment!^!. We obtained the best 
fits to data from 4> asym = —0.6°, 0 asym = —2.5°, and 
P) v =—2 8° 

asym ’ 

Once we have fit these parameters at 4.2 K, the only 
parameters free to fit for the data as a function of tem¬ 
perature are Bo and ui 0 t. For each temperature we sim¬ 
ulated conductivity across a broad range of B 0 and 1 /u c t 
and used a least-squares fitting to arrive at the following 
values for the best fits to data: 


T(K) 

Bo(T) 

7“ct 

4.2 

0 f 

2.50 

14 

2.5 

2.65 

32 

6.8 

3.06 

40 

8.1 

3.28 

50 

8.6 

3.52 

70 

7.2 

4.26 


f We assumed Bo = 0 at 4.2 K in order to perform our fits 
for Fermi surface geometry and alignment. 

The error bars shown on Bq and t/o; 0 r in the main 
text are the standard error of those parameters. At each 
temperature, the values of Bo and V-.T that give the 
best fit to data are those for which the sum of squared 
error (SSE) between data and simulation is minimized. 
We can fit the SSE to a functional form in terms of Bq 
and y UcT about that minimum. We use this functional 
form to approximate the Hessian matrix for these two 
parameters, the inverse of which is the covariance matrix, 
C. The standard error for each parameter is then simply 
given by \/ Cii /(N- 2 ), where N is the number of data- 
points we used for the fitting at that temperature (and 
2 is the number of parameters we fit). Which diagonal 
element of C corresponds to each parameter depends on 
how we construct the Hessian matrix. 


Appendix D: Parameter fitting and error bars 


Appendix E: In-plane transport simulations 


We would expect the parameters k mn to be constant 
with temperature, as they describe the Fermi surface ge¬ 
ometry. Therefore, we can fit these parameters from our 
4.2 K data since we do not expect reconstruction and 
magnetic breakdown to occur at this temperature. From 
band structure calculation^! and from previous AMRO 
studied, we expect this material to have no c-axis dis¬ 
persion along the zone diagonals as well as along the lines 
k x = 7 r and k y = i r. In order for this to be realized, it 
must be the case that 1 — k ei/k 2 i + k wi/k 21 = (P3. 

We simulated conductivity for a wide swath of parame¬ 
ter space and used a least-squares fitting to data to arrive 
at the following: k 0 o = 7.34 nm -1 , Aq 0 = —0.25 nnf 1 , 
kel / k 2 i = 0.69 (and therefore k wi/k 21 = — 0.3lj^4 


In addition to calculating cr zz , we can use the 
same methods as detailed above to calculate the 
in-plane components of the conductivity tensor 53 . 
Neglecting the weak interlayer warping of our 
system, we find v x (ip) = -^kp (ip) cos(cp — 7 ) and 

v v {tp ) = ^rkp((p) sin(i£ — 7 ). Here 7 is the angle be¬ 
tween vf and a vector pointing radially outward towards 
the Fermi surface, and it is given by 


7 (ip) = tan' 


‘ d 
d(p 


(log k F (tp)) 


(El) 


as described in Ref. 1541 The procedure is then nearly 
identical to that for a zz , though slightly simplified by the 
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FIG. 8. (Color online) Three small Fermi surface pockets are formed when the Fermi surface of Tl2Ba2Cu06+i is reconstructed 
under (7r, n) order. 


fact that the A k z ^ terms are not involved in the in-plane 
calculations. We can calculate in-plane conductivity ex¬ 
actly, whereas we can only calculate a zz up to a constant 
of proportionality since we do not know the value of t±. 

Rather than calculating the in-plane transport terms 
and fitting them to experimental data, we want to see 
what predictions we can make for in-plane transport 
based on our analysis of the interlayer transport. We 
fit the points from Figure [5] in the main text to analyti¬ 
cal functions: a second-order polynomial in temperature 
for Y u) c t, and a function of the form ^e~ C2 ^ T for B 0 , as 
we expect that at higher temperatures Bq must decrease 
due to weakening antiferromagnetic correlations. 

Using these analytical functions of our temperature- 
dependent parameters, we are able to calculate the in¬ 
plane transport of Tl 2 Ba 2 Cu 06 +a at any temperature- 
though such calculations should be interpreted with care 
because we are extrapolating to higher temperatures us¬ 
ing information that comes from 50 K and below. We can 
compare these calculations to data taken from compara¬ 
ble samples by Mackenzie et all®, as shown in Figure [ 9 ] 
Note that the data presented in these figures come from 
a sample with T c of 15 K, the same critical temperature 


as the sample whose AMRO data we have analyzed. 

Our simulations of in-plane transport are qualitatively 
similar to experimental data, though they do not agree 
quantitatively, especially the Hall angle and Hall coeffi¬ 
cient. It is important to note that in the magnetic break¬ 
down model, where b o/b plays an important role, we do 
not have Drude-like resistivity: p xy is not directly pro¬ 
portional to the magnetic field. Therefore, we would have 
to use Hall data taken at 45 T to truly make a meaning¬ 
ful comparison. We cannot simply lower the magnetic 
field strength in our calcuations to match the field at 
which data was taken, as we only have information on 
Bq for a 45 T field. It has been proposed that antifer¬ 
romagnetism in the cuprates is enhanced by an applied 
magnetic fiel d 55 * 56 * and therefore we cannot assume that 
the value of B$ at lower fields matches that at 45 T. 

We show these results not because they definitively 
support or contradict the magnetic breakdown model, 
but merely in the spirit of sharing the results of our ex¬ 
plorations. Given that we do not know the dependence of 
B 0 on B 1 it seems unlikely that such in-plane calculations 
can yield strong evidence for or against the suggested 
model. 
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unspecified field). 
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